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ABSTRACT 

We report on the results from two independent but complementary statistical analyses of the WMAP 
first-year data, based on the power spectrum and A-point correlation functions. We focus on large 
and intermediate scales (larger than about 3°) and compare the observed data against Monte Carlo 
ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale 
fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum 
analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high 
at the 95%-99% level (depending on the particular multipole range included). The axis of maximum 
asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but 
tends to lie close to the ecliptic axis. In the iV-point correlation function analysis we focus on the 
northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the 
large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable 
with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry 
is found in the COBE-DMR map, and the axis of maximum asymmetry is close to the one found in 
the WMAP data. 

Subject headings: cosmic microwave background — cosmology: observations — methods: statistical 



1. INTRODUCTION 

In recent months much interest has focused on the 
question of the lack of large-scale power in the WMAP 3 
(Bennett et al. 2003a) first-year data sets (e.g., Efs- 
tathiou 2003a, b; Efstathiou 2004; Spergel et al. 2003; 
Tegmark, de Oliveira-Costa, & Hamilton 2003). Evi- 
dence for non-Gaussian features in the sky maps has also 
been claimed (Coles et al. 2004; Naselsky et al. 2004; 
Vielva et al. 2004; Copi, Huterer, & Starkman 2004). In 
this paper we present evidence of a related effect, namely 
that the large-scale power in WMAP is unevenly dis- 
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tributed on the sky. This effect can be seen on angular 
scales down to 3-5°, or £ < 40. Similar results based 
on the genus statistic have also been reported by Park 
(2004). 

2. DATA AND SIMULATIONS 

All results in this paper are derived from the WMAP 
first year Q-, V- and W-band sky maps, corrected for 
foreground emission using the template method of Ben- 
nett et al. (2003b). The power spectrum analysis in- 
cludes the co-added V- and W-band information, while 
the A-point correlation function analysis utilizes the Q-, 
V- and W-band data both separately and in a co-added 
form. The maps are combined using inverse- variance 
noise weights as described by Hinshaw et al. (2003b). 

The significance of our analyses are evaluated by means 
of Monte Carlo simulations: we generate an ensemble 
of stochastic realizations of the CMB sky to which the 
observed data are compared through the power spec- 
trum and A-point statistics. The CMB simulations are 
based on the WMAP running-index spectrum. From 
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each CMB realization, we generate simulated skies cor- 
responding to the eight WMAP Q-, V- and W-band sky 
maps by convolution with the appropriate beam transfer 
functions (Page et al. 2003) and then adding a corre- 
sponding noise realization, based on the specific noise 
properties of the channel in question. 

3. ANALYSIS METHODS 

In this paper we analyze the WMAP data by means 
of both the angular power spectrum (Hansen, Gorski, & 
Hivon 2002) and its real-space complements, the Appoint 
correlation functions (e.g., Eriksen, Banday, & Gorski 
2002). While we do not aim at presenting the full sta- 
tistical machinery behind these methods here, instead 
referring the interested reader to the above papers for 
details, we do provide the specifics appropriate for the 
WMAP assessment. 

3.1. The angular power spectrum 

The angular power spectrum can be determined using 
a maximum likelihood inference from a local estimate on 
a given patch of the sky, assuming that the coupling be- 
tween harmonic modes due to the incomplete sky cover- 
age can be evaluated. Hansen et al. (2002) have demon- 
strated that for an axisymmetric patch, this correlation 
matrix can be computed analytically. In this paper, we 
utilize a Monte Carlo approach to calculate the matrix 
on a general nonsymmetric sky patch. 

Using this likelihood technique, we have estimated the 
power spectrum on 164 slightly overlapping discs with 
radius 9? 5 uniformly distributed on the part of the sphere 
outside of the WMAP Kp2 sky cut. In particular, we 
have compared the disc spectra in the northern Galactic 
hemisphere to those in the southern sky for the multipole 
bin t = 2 - 63. 

We extend the method further by estimating the power 
spectra directly for hemispheres, thus allowing us to ob- 
tain a higher multipole resolution on large scales (\ow-£) . 
In order to avoid the introduction of bias through con- 
sideration of hemispheres defined only with respect to 
the Galactic coordinate system, we estimated the spec- 
tra for a number of hemispheres defined by orienting the 
north-south axis in 82 different directions on the sphere. 
Specifically, we applied the following procedure both for 
the WMAP co-added V- and W-bands, plus each of 2048 
simulations. 

• The power spectrum was estimated in bins Cb of 
width 3, where the multipole bin b e {CimCax} 
is defined by £ m - m = 2 + 36 and £ max = 4+36. 
The spectrum was successively estimated on the 
northern and southern hemispheres with the axis 
defining north and south oriented in 82 different 
directions. 

• The estimated bins Cb in a given multipole range 
(■ = 4nin — W were summed for each of the 
164 hemispheres giving C h omis P hcrc = Y.b C b- Tnc 
ratio between the northern and southern spec- 
tra for each of the 82 possible orientations r = 
max(CV/Cs, Cs/C/v) was determined. 

• The axis with the maximum asymmetry ratio, 
r max , was recorded. 



Subsequently, histograms of the asymmetry ratio were 
used to quantify the level of agreement between the sim- 
ulations and the real sky measurements. The results are 
presented in H4.ll 

3.2. N -point correlation functions 

An Appoint correlation function is determined from 
the average product of A^ temperatures, as measured in a 
fixed relative orientation defined by an Appoint polygon 
on the sky. Algorithmically speaking, these functions 
are estimated as simple averages over all sets of pixels 
fulfilling the geometric requirements set by the polygon. 

In this paper, the correlation function analyses have 
been designed to probe two different ranges of scales. 

The very largest scales are analyzed by measuring the 
three- and four-point correlation functions from the full 
sky in 1° bins, for a total of 120 (equilateral and pseudo- 
collapsed) 4 three-point configurations and 235 (rhom- 
bic, 2+2-point and 1+3-point; see Eriksen et al. 2002) 
four-point configurations. To facilitate this analysis, the 
large-scale A-point correlation functions are computed 
from low-resolution maps, which are constructed by first 
de-convolving the initial WMAP beam transfer function 
then convolving with a Gaussian beam of FWHM=140'. 
The maps, initially at a HEALPix 5 resolution of A^c = 
512, are then reconstructed at A^ide = 64, applying suit- 
able corrections for the pixel window functions at the two 
resolutions. A best-fit monopole, dipole and quadrupole 
are subtracted from the sky maps before the Appoint 
correlation functions are determined on that part of the 
sky defined by the WMAP KpO mask. However, since 
our processing smooths the maps with a relatively wide 
Gaussian beam, residual, inadequately corrected fore- 
ground structure in the Galactic plane could also be 
smoothed beyond the boundaries of the standard KpO 
mask. For this reason, we extend the mask accordingly in 
all directions before proceeding. This extended mask no 
longer includes the point source exclusion regions since 
the smoothing process minimizes their contribution to 
our analysis. 

The motivation for removing the quadrupole is sim- 
ply that the observed quadrupole amplitude is appar- 
ently low in the WMAP data, and whether this is a real 
feature or simply the outcome of inadequate foreground 
removal remains unresolved (see, e.g., Efstathiou 2003b; 
Efstathiou 2004). Since correlation functions are intrin- 
sically more sensitive to low-^ multipoles and it is dif- 
ficult to include the uncertainties related to foreground 
removal, we choose to eliminate any concern that this 
effect could dominate other interesting results. 

The two-point correlation function is also evaluated 
but it has proved difficult to define an asymmetry pa- 
rameter in this case, and the results are therefore omit- 
ted in this paper. Nevertheless, the four-point functions 
are compared against the value inferred by the two-point 
function of the same sky, as described by Eriksen et al. 
(2002). 

Intermediate scales between 3 and 5° are probed by 
partitioning the sphere into 81 discs, each of 10° radius, 

4 A "pseudo-collapsed" triangle is in this paper defined as an 
isosceles triangle with a baseline between 1 and 2°; see Gaztanaga 
& Wagg (2003). 

5 http:/ /www. eso.org/science/healpix/ 



3 




Fig. 1. — Top panel: The low-resolution Q-, V- and W-band 
co-added WMAP map, to which the extended KpO mask has been 
applied. 

Middle panel: Results from the local power spectrum analysis. 
The color of the large discs indicate the ratio between the I = 2—63 
power spectrum bin of the northern and southern hemispheres, as 
determined in a reference frame where the north pole pierces the 
center of the disc; light red/yellow indicates a low ratio, dark red 
a high ratio. The medium sized dots indicate the absolute value 
of the power spectrum estimated on a 9? 5 disc in the I = 2—63 bin 
(i.e., not on the full hemisphere), and a dark blue dot means that 
this value lies below the lower 80% confidence limit, while a green 
dot means it lies above the upper 80% limit. Finally, the ecliptic 
poles are marked by the small dark-blue spots. The figure should 
not be interpreted as implying that data close to the galactic plane 
have been used in the analysis: the WMAP Kp2 mask has been 
applied for all power spectrum computations. 

Bottom panel: The results from the intermediate scale three-point 
analysis. Blue corresponds to low x 2 > which again corresponds 
to small fluctuations in the three-point function, while red corre- 
sponds to high x 2 - The full KpO mask has been applied to the 
data. 



and computing the three-point correlation function for 
each disc. The full KpO mask is applied to the sky map, 
including the regions of exclusion related to known point 
sources. In order to reduce disc-disc correlations we first 
apply a high-pass filter to the maps, removing all multi- 
poles with £ = 0, . . . , 18. On each of these discs, we com- 
pute the three-point correlation functions for 460 isosce- 
les configurations smaller than 5°. 

The degree of agreement between the simulations and 
the observations are quantified in terms of a standard 
covariance matrix \ 2 statistic. Such a statistic is, in 
principle, only appropriate if the data under considera- 
tion follow a joint Gaussian distribution. Usually, it also 
works quite well for mildly non-Gaussian distributions, 
and in particular symmetric ones, but for strongly asym- 
metric distributions it is likely to yield biased results. It 
is easily seen that the distributions for the even-ordered 
correlation functions for a given geometrical configura- 
tion are in general strongly asymmetric. We therefore 
transform the data of each configuration into a Gaus- 
sian distribution by means of the empirical distribution 
function, before performing the x 2 analysis. The trans- 
formation is defined as follows 

Rank of observed map 1 f s i >i , 

; ; 7- — = —f= / e ~' dt - (!) 

Total number of maps +1 v27r j_oo 

The left-hand side is the fraction of simulations with a 
lower correlation function value than the map under con- 
sideration (i.e., it approximates the true, but unknown, 
cumulative distribution function), and the right-hand 
side yields the corresponding value, s, measured relative 
to a standard normal distribution. 

4. RESULTS 

4.1. Power spectrum 

We have computed a local power spectrum estimate 
for 164 slightly overlapping discs with radius 9? 5, uni- 
formly distributed on the part of the sphere outside of 
the WMAP Kp2 sky cut and compared these to spec- 
tra derived from an ensemble of 6144 simulated maps. 
Concentrating on the lowest multipole bin £ = 2 — 63 
we found that the amplitudes for discs in the northern 
Galactic hemisphere were generally lower in the WMAP 
data than in the simulated maps. Conversely, we found 
that the discs in the southern Galactic hemisphere were 
of generally higher amplitude than in the simulations. By 
considering the ratio of the mean of the spectra in the 
northern hemisphere to that in the southern hemisphere, 
we found that only 0.5% of the simulations have a ratio 
as low as the WMAP data. This is the first evidence of 
a large scale absence of power in one hemisphere of the 
WMAP data. 

In order to pursue this effect further, we have com- 
puted the ratio of the power spectrum amplitudes de- 
termined for the northern and southern hemispheres as 
defined in a particular coordinate system, and for a se- 
lection of multipole ranges. The results are reported in 
Tabled I n particular, we consider this ratio after first 
determining that coordinate frame that maximizes its 
value. The ratio for the WMAP data is larger than 
~ 99% of the maximum asymmetry values determined 
from the simulated maps. We have also tabulated the 
orientation, in Galactic coordinates, of the north pole 
for the data-preferred reference frame. Note that whereas 
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Table 1. Power spectrum asymmetry ratio results 









"eel 






(8, </>) m ax 


WMAP, £ = 


2 


- 40 


0.000 


0.000 


0.003 


(0,0) 


WMAP, I = 


5 


- 40 


0.008 


0.000 


0.007 


(80,57) 


WMAP, i = 


8 


- 40 


0.001 


0.000 


0.047 


(80,57) 


WMAP, I = 


20 


- 40 


0.000 


0.000 


0.009 


(80,57) 


WMAP, I = 


2 


- 19 


0.048 


0.000 


0.002 


(0,0) 


DMR, £ = 


2 


- 19 


0.037 


0.001 


0.131 


(80,95) 



Note. — A summary of results for the ratio of power spec- 
trum amplitudes between the northern and southern hemi- 
spheres over various ^-ranges defined over the WMAP co- 
added V- and W-band data after applying the Kp2 mask. 
First column: The ratio as computed in the ecliptic coordi- 
nate frame. The numbers indicate the fraction of simulations 
with a higher asymmetry ratio for the ecliptic axis than the 
data. 

Second column: The ratio as computed in a coordinate frame 
selected such that the observed asymmetry is maximized for 
the observational data alone. The data values are compared 
against values from the simulations that have the preferred 
axis imposed on them by the data. The numbers indicate the 
fraction of simulations with a higher asymmetry ratio for this 
axis than the data. 

Third column: The ratio as computed in a coordinate frame 
selected such that the observed asymmetry is maximized. The 
numbers indicate the fraction of simulations with a higher 
maximum ratio r than that found in the WMAP or COBE- 
DMR data. Note that in this case the data value is compared 
against values derived from the simulations that may or may 
not have the same preferred axis as the data. 
Fourth column: The (0, <f>) direction of the north pole in the 
galactic reference frame for the axis that maximizes the asym- 
metry observed in the data. 



the asymmetry on the lowest multipoles seems to be con- 
centrated about the north Galactic pole, the asymmetry 
in the higher multipoles (5 < I < 40) seems to be highest 
about the axis with the north pole at (9, (f>) = (80°, 57°) 
in Galactic coordinates. 6 Such a result may argue against 
an explanation in terms of residual foreground contami- 
nation, at least for the higher multipole ranges. It is also 
clear that the observed asymmetry is not simply a re- 
flection of the possible low quadrupole and octopole am- 
plitudes found by the WMAP team. The middle panel 
of Figure ^ summarizes these results in a different way. 
Each observed disc on the map represents the statistical 
deviation, as compared to simulations of the observed 
asymmetry ratio when computed in the reference frame 
for which the north pole pierces the center of the disc. 
What is most immediately evident is that there is an 
apparent lack of large-scale power in the vicinity of the 
north ecliptic pole. 

Figure |2 compares the nearly full-sky power spectrum 
computed by the WMAP team to the local northern and 
southern hemisphere estimates derived in the reference 
frame that maximizes the asymmetry between them for 
the multipole range £ — 5 — 40. This figure also shows the 
best-fit running-index spectrum. We see that the north- 
ern spectrum is systematically lower than the southern 
spectrum over almost the entire multipole range. 

6 Here, 9 and <f> are measured in the HEALPix convention, thus 
corresponding to co-latitude and longitude. 
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Fig. 2. — Power spectra computed from the co-added V- and W- 
band WMAP data. The solid line (histogram) indicates the theo- 
retical best-fit WMAP running index power spectrum. The dashed 
line shows the estimated power spectrum obtained by the WMAP 
team for the Kp2 mask. The black crosses and gray dots represent 
our estimates of the power spectra on the northern and southern 
hemispheres, respectively. Here, north and south are defined with 
respect to the axis that maximizes the asymmetry in the WMAP 
data for the corresponding hemispheres, such that the north pole 
is located at (0,4>) = (80°, 57°). The gray bands indicate the 1 
and 2cr confidence regions, as computed from the ensemble of 2048 
Monte Carlo simulations. Formally, these error bounds differ be- 
tween the hemispheres, but in practice, the difference is small and 
only the values from the northern hemisphere are shown. 



As a useful cross check that should help to mitigate 
against systematic effects as the cause of the observed 
structure, the same hemisphere exercise has been per- 
formed for the co-added 53+90 GHz COBE-DMR map, 
for which we consider multipoles in the range t — 2 — 19 
where the signal is dominant. We find that the DMR 
axis of maximum asymmetry lies close to that for the 
WMAP data. However, for DMR the significance of the 
result is lower at about 87% confidence. Nevertheless, 
given the noisier nature of the data, we consider that 
this is supportive of the asymmetry result. 

4.2. N -point correlation functions 

In Figure|3 the pseudo-collapsed three-point and 1+3 
four-point functions are shown as computed for the 
northern and southern ecliptic hemispheres. While the 
expectation values of the two functions are very differ- 
ent, the observed behavior of the two functions is in- 
ternally consistent: the northern hemisphere correlation 
functions are strikingly featureless (the three-point func- 
tion lies very close to zero, and the four-point function 
drops off very quickly), while the southern hemisphere 
functions show relatively strong fluctuations. 

In order to quantify these statements, we use the full 
covariance matrix x 2 statistic including all bin-to-bin 
correlations. The results from these computations are 
shown in Table [21 The first two rows for each mask (KpO 
and Kp0+|&| > 30deg) indicate the frequency of simula- 
tions with a lower x 2 value than the WMAP data, and 
the third row shows the frequency of simulations with a 
smaller Xn 0rt h/Xsouth- This latter statistic merits some 
explanation. The x 2 statistic in itself measures the over- 
all consistency of an observed function with a predefined 
model relative to the standard deviation of the model. 
Thus, for a function with vanishing mean (such as the 
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Table 2. TV-point correlation function x 2 results 



Mask Q V W Co-added 



Three-point function 



Northern 


KpO 


0.047 


0.023 


0.014 


0.034 


Southern 


KpO 


0.831 


0.861 


0.871 


0.840 


Ratio of x 2 ' s 


KpO 


0.031 


0.014 


0.012 


0.023 


Northern 


\b\ > 30° 


0.046 


0.033 


0.038 


0.041 


Southern 


6 > 30° 


0.792 


0.820 


0.829 


0.793 


Ratio of x 2 ' s 


6 > 30° 


0.040 


0.031 


0.031 


0.039 


Four-point function 


Northern 


KpO 


0.070 


0.054 


0.050 


0.058 


Southern 


KpO 


0.852 


0.873 


0.846 


0.857 


Ratio of x 2 ' s 


KpO 


0.030 


0.021 


0.025 


0.025 


Northern 


|6| > 30° 


0.020 


0.010 


0.012 


0.014 


Southern 


6 > 30° 


0.927 


0.942 


0.931 


0.932 


Ratio of x 2 ' s 


|6 > 30° 


0.004 


0.001 


0.002 


0.002 



Note. — Results from \ 2 tests of the large-scale corre- 
lation functions computed from the northern and southern 
ecliptic hemispheres and from the ratio of northern to south- 
ern x 2 f° r the given bands. Reported here are the frequen- 
cies of simulated realizations with a x 2 value (x 2 ratio) lower 
(smaller) than for the WMAP map. The |6| > 30° constraint 
is enforced in addition to the KpO mask. 



three-point and the reduced four-point functions) the x 2 
statistic is a monotonic function of the overall fluctuation 
level. Therefore, we choose the x 2 itself as a measure of 
the fluctuation amplitudes, and quantify the degree of 
asymmetry between the two hemispheres by a parame- 
ter that is then the ratio of the two x 2 values. 

The conclusions to be drawn from Table El are as fol- 
lows. The northern ecliptic hemisphere shows a consis- 
tently low level of fluctuations, independent of frequency 
and Galactic cut, but inside the allowed range set by the 
model sky. The southern hemisphere shows a slightly 
high level of fluctuations but is again acceptable given 
the model. However, the two are only marginally consis- 
tent internally: as measured by the three-point function 
only ~ 2% of the simulations have a smaller ratio of x 2 
with the KpO mask, and ~ 4% with the additional con- 
straint that |6| > 30°. Somewhat stronger numbers are 
found by the four-point function. Again, about 2% of 
the simulations have a smaller ratio of x 2 m the KpO 
region, but this time only 0.1%-0.4% have a smaller ra- 
tio in the extended region. The latter values tend to 
refute the possibility that residual foregrounds could be 
responsible for the observed asymmetry. Note also that 
the strongest results are found in the bands least affected 
by foregrounds, V and W. 

The top panel of Figure shows the smoothed, co- 
added WMAP map, and in this figure it is actually 
possible to see the asymmetry by eye; while there are 
conspicuous large-scale structures visible on the south- 
ern hemisphere, the northern hemisphere is surprisingly 
featureless. Therefore, in order to localize these features 
further, we have computed the three-point correlation 
function in 460 isosceles configurations smaller than 5° 
on a set of 81 discs, after subtraction of all power on 
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Fig. 3. — Pseudo-collapsed three-point and 1+3 four-point corre- 
lation functions (dotted line) as computed from the low-resolution 
co-added Q-, V- and W-band WMAP data for the KpO Galactic 
mask after imposing an additional constraint that |fj| > 30°. The 
gray bands indicate 1, 2, and 3<r confidence regions as computed 
from a set of 5000 Monte Carlo simulations. The solid line shows 
the median. 



scales £ < 18. The result is shown in the lower panel 
of Figure to aid comparison with the power spectrum 
results shown in the middle panel. An asymmetry is 
also clearly observed on these scales: large connected re- 
gions on the northern hemisphere have a low x 2 ; whereas 
the southern hemisphere regions (especially in the east- 
ern quadrant) demonstrate higher values. The possible 
alignment of structure with the ecliptic poles is in very 
good agreement with both the iV-point results derived 
from the ecliptic hemispheres and the power spectrum 
analysis over the range t = 20 — 40. 

5. CONCLUSIONS 

In this paper we have estimated power spectra and N- 
point correlation functions from the WMAP first-year 
Q-, V- and W-band data, and compared these against 
Monte Carlo ensembles based on Gaussian realizations 
of the WMAP best-fit running-index power spectrum. 
Subsequently, we have found evidence for an anisotropic 
distribution of large-scale power, as implied by the fact 
that the ratio between the northern and southern ecliptic 
hemisphere power spectra is low at the 3<r level. When 
the WMAP data are observed in a reference frame that 
maximizes the asymmetry between hemispheres, the ob- 
served ratio determined from the power spectrum over 
the range I — 2 — 40 is higher than that computed for at 
least 99.7% of the simulations. This is the case even when 
the asymmetry ratio for each simulation is computed in 
its own specific reference frame of maximum asymmetry. 
The significance of the result is therefore independent of 
the reference-frame. The asymmetry was also confirmed 
by an iV-point correlation function analysis, in which the 
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ratio of \ 2 values computed from the ecliptic hemisphere 
correlation functions is low, typically at the 98% level, 
and in some cases, even as low as 0.1%-0.4%. 

The origin of the asymmetry is most simply expressed 
as a remarkable absence of large-scale power in the vicin- 
ity of the north ecliptic pole. Nevertheless, it is also a 
possibility that fundamentally non-Gaussian signals con- 
tribute. Indeed, the three-point correlation function in 
the northern ecliptic hemisphere exhibits what could be 
described as "super-Gaussianity." Such a property had 
previously been determined for the COBE-DMR data 
with a bispectrum statistic (Magueijo 2000). Whether 
this is the case for WMAP remains to be seen, but such 
an unexpected lack of deviation from the zero expecta- 
tion could be interpreted either way. 

It is certain that one criticism of this work could be 
that, given N statistical tests as applied to a specific data 
set, some may prove to be anomalous simply by chance. 
The specious nature of this argument should be clear, 
since the assertion only applies for tests concerning the 
same non-Gaussian degree of freedom, whereas we have 
presented (at least) three. Moreover, it is certainly ren- 
dered moot by the increasing number of recent support- 
ing results in the literature. In particular, Park (2004) 
has independently discovered the north-south asymme- 
try through the use of genus statistics. In addition, the 
low quadrupole and octopole amplitudes already noted 
in the COBE-DMR data have been confirmed by the 
WMAP team itself, but more interestingly, de Oliveira- 
Costa et al. (2004) have noted a conspicuous alignment of 
these features. Finally, several groups have made claims 
for the detection of non-Gaussianity with a variety of 
tests based on phase (Naselsky et al. 2004; Coles et al. 
2004), wavelets (Vielva et al. 2004), and multipole vec- 
tors (Copi et al. 2004). 

That the WMAP data does indeed contain unusual fea- 
tures on large scales seems inescapable. The main issue 
now is to elucidate possible explanations; three possible 
candidates are that these effects are due to systematics, 
foregrounds and exotic physics. 
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It would appear that the WMAP data are remark- 
ably free of systematic artifacts (Hinshaw et al. 2003a). 
Moreover, the results presented in this paper seem un- 
likely to be compromised by systematics since we find 
supporting evidence from similar features in the COBE- 
DMR maps, which are susceptible to different parasitic 
signals. Nevertheless, either previously unknown signals 
could be postulated, or one might accept the possibility 
that the projection of the known systematic signals onto 
these new statistical tests is more complex than expected 
(cf. the COBE-DMR bispectrum situation as discussed 
by Banday, Zaroubi, & Gorski 2000). 

Residual foreground contamination also seems unlikely 
In particular, we note that the t = 20 — 40 power spec- 
trum results find an axis of maximum asymmetry such 
that the north pole of the coordinate system lies away 
from the Galactic plane and the axis separating the hemi- 
spheres is highly inclined to it. More importantly, the 
observed asymmetry is remarkably stable with respect 
to frequency and sky coverage. 

The most intriguing possibility remains therefore that 
the recent findings require some fundamentally new 
physics on large scales in the universe. 
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